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ABSTRACT 


The separated flow past a flat plate at an angle of attack 
has been modelled using the vortex-cloud method The locations of 
the nascent vortices have been determined by applying a simplified 
Kutca condition at the edges of the plate The strengths of the 
nascent vortices have been obtained by taking the average 
shear-layer velocity as determined from the velocities of 
correction of the first four vortices near the plate This 
per-its a feedback from the wake to the separation process The 
calculated flow features are compared with the reported numerical 
results and with experiments 
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CHAPTER 1 


INTRODUCTION 


1 . 1 Discrete Vortex Methods 

The establishment by L Prandtl in 190? that fluid flows at 
large values of Reynolds number can be modelled as essentially 
inviscid flows with thin region of viscous activities near the 
solid boundaries could not be used with profit in flow about bluff 
bodies, where the separation of these thin regions of viscous flow 
(termed as boundary layer) from the solid surfaces disturbed the 
flows rather drastically so that the inviscid methods could not be 
of much use 

The vortex methods for treating such flows evolved slowly and 
built upon the concepts introduced by various researchers 

Th von Karman had shown in 1911 that the wake of a bluff 
body consisted of a series of vortices shed alternately from the 
two separation points Helmholtz has shown earlier in 1938 that 
in an inviscid flow, vortex lines are material 1 1 nes , thati^they are 
continually composed of the same fluid particles and that flows 
with vorticity can be modelled with 'line vortices' of appropriate 
circulation and infinitesimal cross-section It is this 
realization which led to discretization of compact regions of 
vorticity into a collection of vortices embedded in an otherwise 
potential flow 

Rosenhead( 1 931 > appears to be the first person to use this 
approach to study the Ke 1 v i n-He Imho 1 tz instability of vortex 
sheets He discretized a vortex sheet into a number of discrete 
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vortices and allowed these vortices to convect under mutual 
influence with an initial sinusoidal perturbation to the <x,y) 
coordinates of the vortices After a number of time steps vortex 
roll-ups develop along the sheet at the crest and trough of each 
sine-wave culminating finally in a single row of vortex bunches 
(termed as clouds), very similar to half a Karman vortex street 
There is a remarkable and dramatic se 1 f-convect i ve activity 
following a small perturbation 

Rosenhead's success with convection of discrete vortices to 
model vortex sheet had a very seductive influence on later 
researchers Abernathy and Kronauer( 1962) extended this method to 
modelling the alternating vortex shedding in the wake of the 
circular cylinder, a problem of much fascination to researchers 
since von Karman first analyzed it in 1911 Apart from the 

initial di sturbance(F igure 1) the motion is self-induced entirely 
by the convective processes of all the interacting discrete 
vortices Abernathy and Kronauer concluded from this study that 
the role of a bluff body is not very important in the formation of 
its detailed vortex-street wake, excluding its function in the 
generation of the two separate vortex shear layers responsible for 
feeding vortex into the wake 

Much of the recent work on discrete vortex modelling has been 
reviewed by Clements and Maull<1975) Saffman and BakerC1979) and 
SarpkayaC 1989) 

The problem of wake behind bluff bodies has been studied by 
many researchers Almost all such studies involve some variation 
in technique, and all of them rely on some ad-hoc assumption or 


the other 



A/o = 0 281 


Figure 1 Development of von Karman vortex street fro 

sinusoidally perturbed vortex sheet pair (Abernathy and Kronauer 


1962) 
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1.2 Vort-ex -Cloud Approach 


The ’classical’ approach to the problem has been termed as 
the vortex cloud approach for the fact that the discretized vortex 
sheet IS now reoresented by a cloud of discretized point <in 2 
d imens 1 ons ) vor t 1 ces The natural scheme of calculation in such a 

model follows Lagrangian approach wherein discrete vortices which 
are natural points following Helmholtz are convected downstream 
under the combined influence of the inviscid flow and the mutually 
induced velocity of the various vortices constituting the cloud 
A march 1 ng- 1 n-t ime scheme is used The calculations have a 

dramatic advantage in that they are gr i d- i ndependent The 

inviscid flow can be modelled using conformal transformation, and 
so explicit boundary conditions need to be used once the image 
systems are employed to preserve the streamline shapes at the 
body 


This basic scheme has been taken through a lot of variations 
in respect to the manner in which vorticity is introduced 
continually in flow - in the location of introduction of new 
vortices (termed as nascent vortices) and the strength of nascent 
vortices, in the scheme for keeping track of the vortices and in 
methods of cancellation and coalescence of vortices 

C 1 ements < 1 973 ) and Kiya and Arie(1977 1980) prefer the use of 
fixed nascent vortex positions, where new vortices appear 
periodically at fixed location close to (but not at) the 
separation point The rate of vorticity generation is obtained 
from the velocity at these points 


dr 

dt 


= i 

2 • 


( 11 ) 

This formulation is sought to be justified from the pioneering 
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work of Page and Johansen< 1 928 ) on structure of vortex sheets 
Kiya and Arie(1977) allowed the variations of these fixed points 
between subjective limits and compared the calculated results with 
those obtained experimentally to obtain a proper location of 
separation point The effect of the position of the nascent 
vortices on the velocity distribution in the neighborhood of the 
separation point was not examined at all 

SarpkayaC 1 989 ) refers to a work of Kuwahara (Numerical study 
of flow past an inclined flat plate by an inviscid model J Phys 
Soc , Japan, vol 35, 1545(1973)) wherein the strengths of the 
nascent point vortices (at fixed locations) are determined from 
Kutta condition Kuwahara is reported to have found that the 
amplitude of the normal force coefficient calculated in this 
manner depends very strongly upon the assumed position of the two 
fixed points where the nascent vortices appear In fact, a small 
change in the position of nascent vortices is reported to cause an 
eight-fold change in the normal force coefficient 

Based on these considerations, Sarpkaya( 1 975 ) introduced a 
procedure wherein he used variable nascent vortex positions and 


determined the rate of shedding of vorticity from the relation 


ar 

at 



( 12 ) 


where U is the velocity in the shear layer calculated using the 

eh 

average of the convective velocities of the first four vortices in 


each shear layer The location of the appearance of the nascent 


vortices is determined so as to satisfy the Kutta condition at the 


edges of the body This method appears to provide a satisfactory 


feedback mechanism from wake fluctuations to the rate at which 


vorticity IS shed, a shortcoming of almost all the other methods. 
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Sarpkaya working 

with 

others has 

extended 

his 

method 

to various 

other geometries 

like 

cambered 

porous 

plate 

wh 1 ch 

mode 1 

f 1 ow 

through parachute 

canopies rectangular 

pr 1 sms 

and 

even 

about 


oscillating plates(1990. 1986. 1979. 1982) 

The survey of literature of the vortex cloud method shows 
that the method works but only to a limited extent All 
applications involved the use of a number of adjustable parameters 
to make the scheme work and the following main difficulties are 
observed 

a) Groups of discrete vortices moving under their mutual 
influence inevitably lead to a random distribution (Sarpkaya and 
Shoaff. 1979) 

b) Inviscid flow theory does not permit destruction or 
dissipation of vorticity But in real flows only a fraction (say. 


about 60 %) of 

the 

circulation 

fed into the 

shear 

layer is 

found in 

vortex 

c lusters 

(Mair and 

Maull, 1971) 


Different 

ad 

-hoc 

mod 1 f 1 cat i ons 

are 

incorporated 

in models. 

1 n 

violation 

of 

the 

vortex 

theorems 

in order to 

account for 

the 

circulation 

loss 

brought 

about 

by 

the rear face 

shear layer. 

and 

the viscous 

and 


turbulent diffusion 

c) The Lagrangian approach adopted in calculation requires a 
large number of vortices whose velocities due to mutual influence 
need to be calculated This requires n* velocity calculations for 
every time-step, where n is the number of vortices, with n tending 
to be very large Some or the other scheme has to be devised to 
save time Thus, vortex clouds are coalesced somewhat arbitrarily 
to limit computer time 
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d ) Contrary to practical results, the numerical simulation 
of symmetrically started flows rema'n syrwnetrical for a long time 
till round-off errors pile up and start an asymmetry Thus the 

initial time calculation do not give valid results, and large-time 
results are vitiated by the difficulty mentioned in <a> 

above Consequently only a small window of time, large enough to 
start the asymmetrical behavior but small enough for the random 
distribution not having set in, is available for useful results 
This precludes meaningful studies of say opening characteristics 
of parachute canopies 

1*3 Vortex-in-Cell CVIO or Cloud-in-Cell CCIO Methods 

The calculation of velocities of large number of vortices (as 
many as 100 000) by adding the mutual influence of all of them 
stretches the capacity of even the fastest available computer In 
line with the hybrid Lagrang i an-Eu 1 er i an approach of 
Part 1 c 1 e- 1 n-Ce 1 1 (PIC) or Marker- i n-Ce 1 1 (MIC) (Har low , 1 964) a 
Cloud- 1 n-Ce 1 1 or Vortex- i n-Ce 1 1 method has been developed 
(Birdsall and Fuss(1969) Langdon ( 1 970 ) Chr i st lansen ( 1 973 ) , 
Baker(1979), and Smith and Stansby ( 1 988 ) ) These methods combined 
the best features of both the Lagrangian and the Eulerian 
approaches The Lagrangian vortices move through a fixed Eulerian 
gr 1 d 

The method consists of assigning the vorticity of each vortex 
to the four neighboring grid points in the manner of double 
interpolation Thus the vorticity of the n*"^ vortex in a given 
cell is assigned to the four surrounding mesh points (Figure Z ) 
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according to 

W = r A /A^ I = 1 , 2 . 3 , 4 (13) 

V T\ 1, 

where is the circulation of the vortex and A is the area of the 
cell 



Figure 2. Area weighing scheme for the distribution of vorticity 
on the mesh in c loud- i n-ce 1 1 method 

Once the vorticity of all cells is so allocated, the Poisson 
equat i on 

= (14) 

IS solved for the region to obtain for all the mesh 

, j 

points(i,j) The velocity components are then calculated at each 
of the mesh points from the given values of through a simple 

central differencing scheme 


u 

= {y/ 

y> > / 2h 

< 1 3 

) 

V 

4 

1 

11 

- y> ) / 2h 

V-1,J 

( 1 6 > 



where h is the cell dimension 

Once the velocities at the mesh points are obtained, the 
velocity for each vortex is obtained by double interpolation 


u 

n 

W 

It 

A /A 

V 

< 

1 

7 

> 

V 

n 

It 

M 

< 

A /A 

1 

< 

1 

8 

> 


Thus in this hybrid scheme, the Eulerian velocity field is 
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obtained from which the velocity of each vortex is obtained and 
then each vortex is convected in a Lagrangian manner. 

This method requires less time and enables one to track a 
large number of vortices However, SarpkayaC 1989 ) has pointed out 
few drawbacks of this method It makes an otherwise grid-free 
method once again grid-dependent The method and its time-scale 
behavior are sensitive to the size of the mesh, the surface 
boundary conditions, the number of vortices and the time-step 
The flow features of scale smaller than the grid cannot be 
accurately resolved Chr i st lansen ( 1 973 > on the basis of numerical 

experiments has found that the most significant numerical error in 
this scheme arises from the anisotropic CIC - interpolat on of 
velocities 

The present work uses Discrete Vortex Model to s’mulate 
inviscid two-dimensional vortex shedding behind a flat plate It 
uses the method of variable nascent vortex positions as used by 
Sarpkaya< 1973> Numerical experiments are conducted by taking the 
angle of inclination of the plate to the flow of 60 , 85 and 90 
and using the time steps of 0 03 and 0 04 The following c*^apters 
discuss in detail the mathematical description of the model the 
algorithm used and the results obtained 



CHAPTER 2 


FORMULATION 


2.1 The Problem 


This investigation is concerned with the use of the Discrete 
Vortex Method for the study of the separating flow from an 
inclined flat plate kept in a uniform flow with a free stream 
velocity U as shown in Figure 3 The nomenclature used is shown 
in Figure 4 



Figure 3 The wake behind a flat plate at 45° angle of attack at a 
Reynolds number of 4300 Aluminum flakes suspended in water show 
its characteristic sinuous form ( Reproduced from Annual Review of 
Fluid Mechanics, Vol 13, 1981 ) 



Figure 4 Flat plate at an angle of attack 
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2.2 Mat-hema-tical Descrip t-ion 

2.2.1 Conqslex Pot-ential and t-he Transformation Used 

The approach used consists of transforming the plate plane z 
into a circle plane using Joukowsky transformation This relates 
the to z-plane by the relation 

z 

z = C + - ( Z 1 > 

C 

The complex potential for flow in the circle plane at an 
angle c* as shown in Figure 5 can be written as 


w<C> = U C e 


, , 2 ict 

Dee 


( 2 Z > 


where c is the quarter-chord of the flat plate or the radius of 
the circle in the C“Dlane 



Figxjre 5. Flow in the Ca> physical <plate> plane, and 
(b> transformed (circle) plane 
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This represents a non-separate flow past the circle ( or flat 
plate in the z-plane > To si mu late the vortex sheet of the 
separated flow, we add a senes of discrete point vortices in the 
flow Let us use the following nomenclature Let F represent 

kp 

the strength of the counterclockwise k^^ vortex shed from the 

trailing edge ( x = f 2c > of the flat plate in the z-plane and 

^kq strength of the clockwise vortex shed from the leading 

edge ( x = - 2c ) Uie differentiate between the two ser es for 

ease in coding including the fact that with this approach we can 

take both F^^ and F^^ to be positive throughout 

The complex potential for a p-vortex, therefore is 

F. 


kp 

277 


in ( C - 


and that for a q-vortex is 

F 


- i 


kg 

2tt 


In ( C - C,, > 


But these potentials can not be added directly to that of 
flow past a circle, since it would change the shape of the 
stagnation streamline which would no longer satisfy the boundary 
condition of the flow 

For this purpose, we use the circle theorem For a vortex of 
strength T /Zn at a point C = ? + i i}' outside the circle lCl=c, 
the circle theorem gives the complex potential as 


-<C) = . ^ in < C - C' > - . ^ in ( C - ^ in C 

{ 2. ■5 > 

where the first term on the right represents the potential of the 

i 

vortex situated at C and the second term that of the inage of 
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this vortex located at the inverse point and having a strength 
- r/Zn (Ml Ine-Thomson , 1969) 

The last term represents the image at the centre But in the 
present case, there are no images at the centre of the cylinder 
because the vortices have been shed from the cylinder and leave 
circulation opposite to their own on the body This last term is 
equal and opposite to the circulation due to vortex at the centre 
and thus cancels it 

Making use of Equation Z 3, the complex potential function w 
in the circle (^ ) plane (see Figure £r> which describes a uniform 
flow U with a doublet at the origin to simulate the cylinder, k+1 
vortices rotating counterclockwise on the right-hand side of the 
wake (called p-vortices), k+1 vortices rotating clockwise on the 
left-hand side of the wake (called q-vortices>, and the images of 
all the p- and q- vortices in the circle may be written as 


w(C> = 


- U ( C + s'" > + 1 ^ in ( C - C > 


Zn 


- 1 


op 


In ( C - r ^ ^ I In ( C - 


Zn 


r k=i Zn 

op 




2 r 

In ( < - I ) - I In ( ^ 

k=i Zn C, 


m r 
- Z 1 - 


kp 


Zn 

r. 




+ I 


r* 2 m 

In ( C - - > - 2: 1 In ( C 

C k = i 

o<i 


Zn 

m r 


Zn 


- 


+ z , JS3. In < c - § ) 


k = i Zn 


kq 


( 24 ) 


in which r and C represent respectively the strength and 

kp kp 
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location of the p-vortex, F, and f, 

kq kq 

location of the q-vortex, and c the radius 

an overbar indicates a complex conjugate 


the strength and 
of the cy 1 i nder 


2* 2# 2 Nascent Vortices and the Kxitta Condition 
r and r located at r and C respectively are called 

op oq op oq 

the ^nascent vortices* or ^Kutta vortices’ These represent the 
first vortices introduced into the flow at each edge of the plate 
The positions of these nascent vortices are chosen so as to 
satisfy the Kutta condition at each edge of the plate In 
accordance with this condition, the flow separates tangentially at 
the edges of the plate Thus, 

dw/d'^ = 0atC==~c <25> 

in the circle plane and dw/dz gives a finite velocity at z = ± Zc 
in the plate plane Inserting Equation 24 in 25, we get two 
equat i ons 


r . 

, , , ^1,01 ^ , <=>p I 

- U ( e - e ) f 1 — — 


- i 


op 


+ s: I 


kp 1 


kp 


f I 


oq 


Zn 

m 

Z 1 

c - 

C 2^ 

op 

1 

Zn 

2 ~ 

k=l 

c - c /c, 

kp 

m 

r 

_ m 

Z 1 

kq 

1 + z 

k=± 

Zn 

c - C,, ■<=* 


2 


op 


oq 


Zn c - c 


kq 


oq 


1 


oq 


= 0 


( 26 ) 
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U ( - e"® > - 1 _i + 1 ^ L 


Zn c + C 


op 


Zn c + /C 


op 


Z i 


kp 1 


+ Z i 


kp 1 


+ 1 


oq 


k=i Zn c + C, k=i Zn c + c /C, 

kp ^ kp 


Zn c + C 


oq 


- I 


oq 


1 


+ Z I 


kq 1 


Z I 


kq 


Zn c + C^/C k=i Zn c + C, k = i Zn c + c^/C, 

oq kq kq 


= 0 


( 2 7 > 


Equations 2 6 and 2 7 are solved to get the positions of the 
two nascent vortices While solving these equations, two 
assumptions are made for simplification 

<i) The vortices T and F are introduced alonq the 

o p oq 

^ -ax IS, 1 e , 


77 = 77 = 0 < Z 8 > 

op oq 


(ii) The effect of F and F on each other at the time of 

op oq 

introduction is negligible 

Thus in Equation Z 6, we neglect the term containing F and 

oq 

Similarly in Z 7, the term containing F 

op 

Simplifying the two equations with the above considerations 
we get the positions as 


and 




2 


c 

C 

i. 

F 

op 

/ Zn 

1 + 

o 

o 

F / 2rr 

op 


2 

c c 

2 

W 

r / 27 t 

oq 


c c 

2 

X » 

r / 2tt 

oq 


< 2 9 ) 


( 2 10 ) 
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whei'e 


c = -ZUsin a 
1 


Z [ 

k = 1 Zn 


c + ? 


kp 




■] + ^ 


kq 


C + ? 


kq 


k = i Ztt 


c < c - ? >■ 

kq 


( 2 11 > 


c = -2Lls 1 n a + Z 
2 

k=i Zn 


kp 


c - ? 


kp 


c < c + ? > 

kp 


] - Z 


kq 


kq 


k=± Zn 


c ( c + ? )• 

kq 


( 2 12 ) 


2.2.3 Velocities in the Circle and Plate Planes 

The advancement of the position of each vortex in the z-plane 
at the end of each time interval requires the transformation of 
the positions of all vortices into the circle plane, the 
calculation of the velocities at each vortex center in the circle 
plane and, finally the reverse transformation of the velocities 
from the circle to the plate plane 

For calculating the velocity components in the circle plane 
of any vortex the say, the complex pote'^tial of that vortex 

is subtracted from the total potential giver by Equation 2 4 and 
evaluating the derivative of the remaining terms at C = Cj^ 

Using the above procedure, the potentials of the p-vortex 

t-h 

and the k q-vortex may be written as 


w, ( C > = w ( C > 

kp 




( 2 13 > 


and 




w < C > + 



in < C - > 


( 2 14 > 
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The velocities, then, in the C'Pla^ne are given by 


<Pk 

? V 


and 


[ / dC ^ ^ 


< 2 15 > 


kp 


? u 


[ dWj^^<C> / dC ^ ^ 


< 2 16 > 


kq 


The velocities in the plate plane are then calculated using 
the relations 


kp kp 

X y 


< -U, + 1 V. ) 

kp- 


kp r-Z 2 

>> '^kp- = 


r r ^ 

kp ^ kp kp 


2 rT 2.2 


( 2 17 > 


- u + 1 V 

kq kq 

X y 


( - U, + I V, > 

kq- 


kq -2 2 

^kq- ^ 


^2 _ ^ 2 
^kq ^ ^kq ^kq 


Zn ,_2 2.2 

"^k,- = ' 


( 2 18 > 


which are derived usina the transformation 


2.3 Method of Calcxilation 
2.3.1 Vorticity 

The rate of shed of vorticity into the wake is given by 

where <5 is the boundary layer thickness This may be approximated 
by 

^ = 1 ( V® - ) ( 2 20 > 

dt 2 12 

where V and V reoresent the velocities at the outer and inner 


edges of the shear layer 
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Page and Johansen <1928), on the basis of experiments with 
various bluff bodies have shown that vorticity is shed from the 
two sides of an asymmetric body at the same rate and that the 
velocity at the outer edge of the sheet is much larger than the 
velocity at the inner edge Thus we can ignore the term in 
the above expression for vorticity flux to simplify the expression 


as 



< 2 21 > 


Now the question comes of selecting the velocity There has 

been varied opinion by investigators regarding the selection of a 

point at which the velocity should be considered for calculating 

the vorticity flux as also mentioned in chapter 1 Cl ements ( 1 973 ) 

and Clements & Maull<1973> selected a fixed point in the flow near 

the separation point and used the velocity U at that point to 

8 

calculate the rate of shed of vorticity using Equation 1 1 

Kiya and Arie(1980) proposed a 'velocity-point scheme' , in 
which a 'velocity-point' is introduced in front of the plate along 
a line perpendicular to the plane of the plate and drawn through 
the edge concerned, which they interpreted as the edge of the 
boundary layer at the separation points The velocity U at that 

3 

point was used to determine the strength of the nascent vortices 

Sarpkaya ( 1 973 , 1 989 ) , however, felt that in these methods no 
interaction is allowed between the shed vortices and the amplitude 
of oscillation of the point or the time of appearance of the 
nascent vortices He used the velocity U m the shear layer for 

3n 

calculating the vorticity flux (Equation 1 2> 

It IS felt that the various schemes are quite arbitrary In 
the present work. Sarpkaya ' s ( 1 975 ) approach has been used as it 
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appears to simulate the mechanism of feedback from wake 
fluctuations to the fluctuations in the rate o- circulation 

U IS calculated by using the average of the transport 

3h 

velocities of the first four vortices in each shear layer Thus 


4 . 1/2 

U ^ i Z ( t v^ > 

% k=i "'p 


4 1/2 

U^=-Z(u, +v, ) 

ah 4. kq kq 

q k=l 


< 2 Z2 ) 

( Z Z7> ) 


Sarpkaya ( 1 975 ) suggests that the use of t^e velocity at the 

position of a single vortex (with k=l say) to estimate the shear 

layer velocity would not be satisfactory parti f because or F^^ 

vortices are too close to the edges of the plate and the radii of 

curvature of their pathlines are relatively small Thus the 

numerical errors involved in the determination of their transport 

velocities could give rise to unacceptable oscillations in the 

vorticity flux The use of the average of the transport 

velocities of a number of vortices in each snear layer would 

therefore give a better representation of the velocity and 

vorticity flux in the shear layers After several preliminary 

calculations, he has concluded that taking U ^ as the average of 

sn 

the velocities of the first four vortices will be quite 
sat i sfactory 


2.3.2 Int.roduct.ion and Convect-ion of Vortices 

The method of introduction, as well as the subsequent 
convection, of the nascent vortices plays an iiwDortant role in the 
evolution of wake by discrete vortex method They are required to 
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be 1 ntroduceci at very small time intervals and convected in a 

manner that they can follow the true streamlines in the regions of 

flow where the streamlines are rapidly curving 

To introduce the first vortex, use is made of the results of 

Page and Johansen ( 1 928 ) who showed that the velocity in a shear 

layer in a steady flow is about 1 4U , where U is the freestream 

velocity Thus from Equation 2 21, we get 

dr ^ ± 2 

dt 2 1 

^-(14 U>^ 

2 

^ ( 2 24 ) 


As the non-dimens 1 onal free stream velocity is equal to 
unity we have 

AT ^ At < 2 25 ) 

Accordinaly, the strenaths F and F of the first vortices 

“ - Ip iq 

are taken equal to 5 At The positions of F and F are found 

±p Iq 

using the Kutta condition Equations 2 6 and 2 7 which are reduced 
with only the first vortices present in the flow to 


F , F , 

,,, -vO( tot. ip 1 , ‘P ^ 

U(e - e ) + i -=-4- - i 
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ip 
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iq 1 

<C + > 

iq 


0 

( 2 26 ) 


0 




< 2 27 ) 
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Here an assumotion is made that these starting vortices are 
introduced along the ?-axis, i e , 


r). 




0 


( 2 28 > 


■ip 'iq 

This assumotion simplifies the equations to a great extent and 
the ^-coordinates o-'' the first p- and q- vortices can be found 
solving the simple algebraic equations 

Now, to calculate the velocities of these -^irst vortices, we 


write the potential for the situation when only F and F are 

ip i<3 


present 

2 F 

r X It r '■“V’Ol , C 'uOl ^ 

w(C ) = -U(Ce + — e > + i 


F 2 

^ ln(C - C > - I ^ InCC - — ) 
2tt ^ ^ ip Zn ^ = 


- I ^ ln<C - C > + 1 InCC - ^ > 

2tt Iq ^TZ r 

( 2 29 ) 


w <C> = w<C> - 1 ^ ln<C - C, > 

ip ^ ^ Zrz ip 


2 

. , - la , c vO(. 

- U (C e + — ® ) - i 

c 


F 2 

ln<C - > 

ip 


' ^ ^ ‘ ^ - f- > 


( 2 ^0 ) 


From this, the velocity of the first p-vortex can be derived as 


-u + 1 V 
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Solving the right-hand side of this equation, the rea’ part 
gives the ^-component and imaginary the 7 j-component of velocity of 

r 

ip 

In the physical plane 


< - u + i V > 
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X y 
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Thus, after these calculations, we know the positions and 


velocities of the first vortices introduced into the flow 
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After introducing T and F 

ip iq 


vortices AF and AF 

Op oq 


we introduce two nascent 
each of strength equal to At in accordance 

with Equation 2 25 

First vortices are then convected in physical plane for a 
time interval of At in five steps of — At using the expressions 


x(t + i At) 

5 

y<t + - At) 

5 


x(t) + - u(t) At 

!5 

y(t) + i v(t) At 

5 


< 2 36 ' 

i Z 37 ) 


After the time step At, we again add AF = AF = At This 

op oq 

procedure is repeated five times till the first vortices are 
convected for 5At The velocities of F and F vortices ar® 

1 o 

computed after every time interval At, and thus the undated 

velocities are used for convection 

At time 5At , the last nascent vortices are renamed as first 

vortices and the first vortices are renamed as second We repeat 

the above procedure of introduction and convection four times in 

order to have four vortices in each shear layer 

Having introduced four vortices of constant strength we can 

now introduce new vortices of strengths calculated using the 

velocities and defined earlier and given by Equations 

p q 

2 22 and 2 23 Thus, 


AF 

1 

At 

op 

2 

oh 



P 

AF 

1 

U^,. At 

oq 

2 

oK 


( 2 38 ) 

< 2 39 ) 


Positions C and C ot the nascent vortices F and F are 

Op oq op oq 

then calculated using the Kutta condition equations. 

We then convect the first vortices in 3 steps of — At and 

5 
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also calculate their velocities prior to each steo 

The velocities of rest of the vortices are then calculated 
and these vortices are convected in a single steo of At 

We repeat the above steps 5 times so as to convect a’ 1 the 
vortices for 5At 

At this stage, the vortices are renamed and new •'ascent 
vortices are introduced into the flow 

This way we evolve the wake behind a flat plate for ar .• time 
period by introducing the vortices accordingly 

2.3.3 Removal and Coalescence of Vort-ices 

The vortices which are very close to the plate induce very 
large velocities and lead to absurd results Sarpkaya< 1 97 ' ) has 
argued that vortices which approach the plate will tend to 
dissipate because of viscous action close to plate and, the-'efore, 
it might be prudent to remove such vortices A scheme has been 
used wherein any vortex that comes within 0 Ic (2 5% of the plate 
length) is removed from calculation 

K 1 ya and Arie(1980) have also discussed this p-'oblem 
According to them, it is probable that individual vortices come 
too close to the rear face of the plate and this causes t''em to 
have unreasonably high velocities along the plate owing to the 
presence of the image vortices within the circle !•' the 
transformed plane He , therefore , also removed the vortices in his 
model whenever they came nearer to the rear face of the plate than 
a specified distance 

It has been well recognized that the orbiting motion o-^ point 
vortices can rapidly affect neighboring vortices and that the 
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vortices of opposite circulation could remove themselves from the 
wake at high speeds when in close proximity This is avoided by 
coalescing the vortices Two vortices having same sign (i e 
positive or negative) are coalesced into a single vortex of the 

strength equal to the sum of the two whenever they come closer to 

each other than a distance equal to 0 Ic In a similar manner 
two vortices of opposite sign also are removed from calculation 
whenever they come closer than 0 Ic to each other The distance 

0 Ic has been selected from the observation of Schaefer and 

Esk 1 naz 1 ( 1 959 ) who reported that a line vortex which has been in 
the wake of a viscous flow with a Reynolds number of about 40 000 

for a period t=8c/U will have a core radius of approximately 

0 07c Thus the distance 0 Ic is of the order of the core radius 

of a decaying vortex 

By the above methods, quite a large number of vortices are 

removed and the number of vortices to be handled decreases saving 


a lot of computer time 



CHAPTER 3 


RESULTS AND DISCUSSION 

The computer code for the present work has been made in 
FORTRAN 77 language and run with the desired input parameters on 
the Convex 220 system 

Numerical experiments were done to evolve the wake of a flat 
plate immersed in a uniform flow Keeping in mind the 

availability of computer time computations were made for only 
three cases as far as the angle of attack is concerned, viz for 
flow normal to the plate, inclined at 85"* to the plate and at 60° 
to the plate To start with, the time-steo At of 0 04 was 
selected for convection of vortices following Sarpkaya< 1 975 > who 
on the basis of a number of numerical experiments with various 
values of At has concluded that this is the optimum time-step 
compromising between the computer time and the accuracy of 
results The computations were done, however with one more value 
of At , 1 e , 0 03 

The results for all the cases studied are summarized below 

3. 1 Case !♦ Plat-e Normal t.o the Flow 

3.1.1 Time-step At = 0.04 

The program was run for t = 180 The results plotted in 

Figure 6 for various intermediate values of time starting from t = 
10 show the step-by-step evolution of the wake From the figure, 
it can be seen that the vortices start shedding symmetrically from 
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both the edges of the plate and curve inwards in a very well 
defined manner This type of symmetric pattern, however, does not 
last for long, and after t = 20, the asymmetry grows and a phase 
difference starts building up between the p- and q- vortex 
clusters At t = 70, the picture clearly resembles the uniformly 
accelerated flow picture around an airfoil at 80° angle of attack 
(Figure 7 At around t = 80, vortex clouds start shedding from 

both the edges of the plate The majority of the vortices 

congregate in clouds and move with the clouds as a whole At the 
later time the distribution of these vortex clusters clearly 
resembles the von karman vortex street The availability of 

computer time put a constraint on running the program further 
The compu*er time taken by the run was 8 hours 

The - ' ow pattern thus clearly exhibits the following major 
features as also noted by other investigators 

i) ”ne wake starts developing symmetrically and it is only 
after some time that the asymmetry sets in Unlike the work of 
Sarpkaya and Shoaf f ( 1 979 ) , who used artificial perturbation method 
to induce asymmetry while modelling the flow over a circular 
cylinder in the present computation , the flow became asymmetric 
by its own It is felt that the rounding errors probably have 
caused slight asymmetry, which once set in, grows with time 

In this context, it is appropriate to report the work of Kiya 
and Arie (1980) who inclined the flat plate by 5“ from the normal 
position to avoid the introduction of artificial asymmetry of flow 
to initiate the periodic vortex shedding 
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Figure 7 Start of asyovnetry of wake 

a) uniformly accelerated flow around an airfoil at 80“ an 
attack (Reproduced from American Scientist, Vol 72, 1984), 

b) wake flow pattern behind the flat plate during the early 
of flow development computed in the present study 
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Figure 8 Vortex arrangements behind a flat plate c* = 50 

(Sarpkaya, 1975) 
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Figure 9 Vortex pattern behind a flat plate 
Arie, 1977) 
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n> Once the asymmetry sets in. the pattern of vo-tex 
shedding is established and we can make an estimate Strouhal 

number 

St = f L/U 

where f is the frequency of vortex shedding 

In this phase of regular vortex shedding. the frequency is 
estimated as about 1 in 25 seconds giving a value of St as 0 16 
This compares favorably with the experimental results of 
Abernathy ( 1962) who obtained St = 0 164 in the range of 40“<af80° 
and the calculated results of Sarpkaya< 1 975 ) for the same range 
and Kiya and Arie(1977) (Figures 8 and 9> who calculated the value 
of St as 0 154 and 0 14 respectively 

3.1.2 Time-step At = 0.03 

The program with At = 003 could be run for t = 130 In this 
case also the flow started symmetrically but became asymmetric at 
a later stage (Figure 10) This time, the periodic vortex shedding 
started at time t = 60, which is lesser than with At = 0 04 This 
IS justified considering the fact that the rounding off errors are 
introduced faster with smaller value of time-step 

The vortex shedding frequency, however, in this case was 1/15 
as against 1/25 for At = 0 04, giving St = 0 26, which is much 

higher than the expected result 

3.2 Plate xnclined at 85° to the Flow 

For <a = 85*, the program was run with At = 0 04 only and for 
a time t = 70 The results for c» = 85 are plotted in Figure 11 
The asymmetry, in this case starts right from the beginning and 
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Figure 10 Evolution of wake behind a flat plate c* = 90 , 
0 05 (contd > <e) t = 60. (f) t = 70. (g) t = 80 
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a = 85“. At = 0 04 (a) t = 10, (b) t = 20. <c) t = 50, (d) t 
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Figure 11 Vortex shedding pattern behind an inclined flat plate 
a = S?* , At = 0 04 (contd > <e) t = 50, (f) t = 60, (g) t = 70 
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vortex clouds shed alternately It was Kowever, difficult to 
find out the Strouhal number in this case from the vortex shecding 
pattern as the frequency of shedding could not be establ shed 
from data output 

3.3 Plat© Inclined at 60“ to the Flow 

The vortices in this case start shedding asymmetrical! as 
expected The picture, however, was noticeably different fro- the 
other cases described earlier Here, the p- vortices which are 
shed fromthe trailing edge at the start of the flow leave the edge 
and convect downstream and for a long time the other o- vc'tex 
cloud IS not formed The q- vortex cluster on the leading edge 
grows in size during this period It is only at t = 170 that the 
second p- vortex cloud starts forming At t = 180 the c ouds 
leave the plate and shed with a phase difference as can be seen 
from Figure 12 

This type of delayed vortex shedding, however, is not 
considered to be due to any error in the code as Kiya and 
Arie(1977> have also reported that the periodic flow patterns 
cannot be predicted from vortex patterns in an early stage of the 
flow development 

As the program could not be run beyond this time, the flow 
features could not be studied fully and the Strouhal number could 
not be estimated 

It was tried to compute the results for this angle of attack 
with a time step of 0 05 also The program, however, could net be 
run for a longer time The results for t = 10 to 40 are plotted 
in Figure 13 
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Figure 12 Vortex pattern behind an inclined flat plate a = 60 , 
At = 0 04 <a) t = 10. <b> t = 20, (c) t = 30 
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Figure IZ Vortex pattern behind an inclined flat plate a 
At = 0 04 (contd ) (d) t = 180, <e) t = 190 
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Figure 13 Vortex pattern behind an inclined flat plate a = 60°, 
At = 0 03 (a) t = 10, <b) t = ZO, (c) t = 30. <d) t = 40 
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On analyzing the data for various cases studied it is felt 
that probably a little more attention has to be paid to study the 
way the nascent vortices are introduced and the Kutta condition 


CHAPTER 4 


CONCLUSIONS 


In the present work the vortex shedding behind a flat plate 
has been numerically studied using the discrete vortex method, in 
which the shear layers which emanate from the leading and trailing 
edges of the plate are represented by an array of discrete 
vortices through the use of the complex potential and the 
conformal transformation between the plate and the circle plane 
Kutta condition has been used to determine the locations of the 
nascent vortices The model which is a simpler implementation 
than the other reported methods reproduces the results that have 
been reported The onset of asymmetry occur earlier in our 
procedure, probably because of the approximation involved in the 
application of Kutta condition 

The Strouhal numbers predicted are around the same values as 
reported earl i er 

This implementation has not been able to calculate 
satisfactorily the flows for small values of angle of attack 
wherein the chaos sets in before recognizable vortex clusters are 
developed This may be because of rather severe dissipation of 
vortices in the rear-face shear layer modelled here 

It IS felt that to save building up of truncation errors, and 
incidentally to save time for calculations, the vortex clusters 
once they are identifiable should be replaced by a coalesced 
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single vortex SarpkayaC 1975) and Kiya and ArieC1977) 

implemented such schemes with benefit The scheme could not 
implemented in the present study due to lack of time 

Another apparently useful approach to tr/ is to combine 
vortex cloud method with the c loud- i n-ce 1 1 method None of 
reported work has handled the separation phenomena 

c 1 oud- 1 n- c e 1 1 method One could use the vortex cloud method 

the formation of the nascent vortices, but onte the vortices 
formed, the convection can possibly be handled by 
c 1 oud- i n-ce 1 1 procedure 
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